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I.  DESCRIPTION  OF  THE  PROJECT. 


A.  Desired  result. 

The  final  goal  is  to  devise  means  to  maintain  continuous  contact  with  the  hypersonic 
vehicle.  When  such  vehicles  were  principally  spacecrafts,  a  blackout  period  of  up  to  two 
minutes  was  acceptable  albeit  undesirable.  But  when  the  vehicles  are  of  military  origin,  it 
is  clear  that  continuous  contact  is  essential  for  both  targeting  and  rapid  abort  reasons. 

B.  Challenge. 

A  vehicle  moving  through  the  stratosphere  (altitudes  40km-50km)  at  hypersonic  velocities 
(8-15  Mach)  is  covered  by  a  plasma  sheath.  Typically,  the  plasma  density  n  can  be  as  high 
as  1018m-3  with  corresponding  plasma  frequency 


(1) 


of  about  9 GHz.  In  (1),  e  is  the  electron  charge  — 1.6  x  10  19C,  eo  =  8.85  x  10  12CV  1m  1 


and  M  is  the  electron  mass  9  x  10  31  kg.  Therefore  the  plasma  is  opaque  to  frequencies  lower 


than  9 GHz.  Direct  communication  through  such  a  plasma  to  and  from  the  vehicle  is  impos¬ 


sible  because  frequencies  /  suitable  for  long  distance  propagation  through  the  atmosphere 


are  usually  much  less.  For  example,  the  standard  frequency  used  for  navigational  satellite 
systems,  including  the  global  positioning  system  (GPS),  are  less  than  2 GHz.  For  the  GPS, 
/  =  1.57542G#*. 

C.  Review  of  solutions. 

It  is  a  challenge  which  has  drawn  many  responses.  They  fall  into  several  categories. 
The  first  ignores  the  presence  of  the  plasma  by  using  signals  with  frequencies  well  above  the 
plasma  frequency.  The  difficulty  with  this  method  is  that  such  signals  are  heavily  attenuated 
in  and  scattered  by  the  atmosphere.  A  second  means,  which  also  ignores  the  plasma,  is  to 
use  low  frequency  signals  in  the  100 MHz  range  where  wavelengths  are  large  compared  to 
the  plasma  sheath  thickness  (typically  of  the  order  of  a  meter).  But  such  solutions  have 
high  cost  and  low  bit  rates  and  are  not  well  supported  by  existing  infrastructure.  A  third 
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category  of  solutions  violates  the  plasma.  One  approach  is  to  remove,  by  vehicle  reshaping, 
for  example,  the  plasma  from  certain  points  on  the  vehicle  at  which  one  might  place  an 
antenna.  Another  is  to  destroy  it  by  electrophilic  injection  or  by  injecting  water  drops.  A 
third  approach  is  to  use  powerful  magnets  to  reshape  the  plasma.  Such  solutions  involve  a 
heavy  cost  in  that  design  features  necessary  for  their  implementation  must  be  built  into  the 
vehicle  a  priori.  Nevertheless  some  are  feasible  and  worthy  of  consideration.  For  example, 
it  is  possible  to  build  an  antenna  into  a  sharp  leading  edge  which  would  protrude  beyond 
the  plasma  and  survive  for  sufficiently  long  (it  would  be  eventually  destroyed  by  ablation) 
to  cover  the  flight  time. 

D.  The  most  attractive  and  perspective  solution. 

The  fourth  category  of  solutions,  and  the  one  to  which  we  are  attracted,  uses  the  prop¬ 
erties  of  the  plasma  itself  to  effect  transmission  in  the  same  way  a  judo  expert  uses  the 
strength  and  motion  of  an  opponent  to  defeat  him.  One  idea  is  to  create  new  modes  of  os¬ 
cillation  and  propagation  by  the  introduction  of  magnetic  fields.  Indeed,  for  strong  enough 
fields,  the  Larmor  frequency  fLarmor  is  sufficiently  large  that  the  window  (fLarmor,  max(/^)) 
for  which  the  plasma  is  opaque  is  small  and  transmission  can  be  achieved  for  frequencies 
below  fLarmor ■  But  the  introduction  of  magnetic  fields  involves  large  additional  weight  and 
new  design  features.  The  second  idea  is  much  more  simple.  Its  aim  is  to  take  advantages  of 
nonlinear  properties  of  plasma  to  render  it  effectively  transparent  to  the  signal.  Communi¬ 
cations  both  to  and  from  the  vehicle  are  feasible  using  basically  the  same  ideas.  We  shall 
first  describe  the  “to  the  vehicle”  case.  Consider  Figure  1  in  which  we  show  schematically 
the  response  of  the  plasma  to  an  incoming  signal  with  low  frequency  u>  from  a  direction 
which  makes  an  angle  (p  with  the  normal  to  the  vehicle.  There  are  two  principal  features  to 
the  response.  First,  there  is  a  reflection  from  the  layer  at  a  point  z  =  zr  where  the  plasma 
frequency  at  the  point  lol(zt)  is  u>  cos  <f>.  However,  the  influence  of  the  signal  is  felt  beyond 
that  point,  namely  at  the  resonant  layer  2  =  0  where  cjl(0)  =  ui.  Langmuir  oscillations  are 
excited  there  which  produce  large  transversal  and  longitudinal  components  of  the  electric 
field.  The  resonant  layer  acts  as  an  antenna.  The  task  is  to  find  a  way  to  connect  the 
antenna  at  the  resonant  layer  at  2  =  0  to  a  receiver  on  board  the  vehicle  at  2  =  R.  There 
are  several  possibilities  which  we  have  outlined  before  [1-3]. 
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FIG.  1:  u>L(zr)  =  uj  cos  (p, u>l{0)  =  id.  If  the  thickness  of  the  plasma  sheath  is  equal  to  L  +  R  =  1  m, 
the  signal  frequency  /  =  2 GHz  and  the  plasma  frequency  ft  ~  9 GHz  then  L  ~  5 cm  and  R  ~ 
95  cm. 

The  most  practical  one,  however,  is  also  the  most  simple  and  first  suggested  without  a 
detailed  numerical  simulation  in  [1].  We  use  an  onboard  source,  which  we  call  the  pump, 
to  generate  electromagnetic  signals  of  sufficiently  high  frequency  uip  (ujp  >  ma xul(z)  +  u>) 
that  they  can  propagate  through  the  plasma.  There  are  several  candidates  for  such  a  source. 
For  example,  available  on  the  open  market  is  a  klystron  amplifier  which  can  generate  3 kW 
of  power  at  frequencies  of  12  —  14 GHz.  These  high  frequency  waves  have  only  to  travel 
distances  of  a  meter  or  less.  They  interact  nonlinearly  with  and  scatter  off  the  signal  wave. 
Not  surprisingly,  the  largest  contribution  to  the  scattered  wave  comes  from  the  nonlinear 
interaction  of  the  pump  wave  with  the  plasma  density  distortion  induced  by  the  incoming 
signal  wave  at  the  resonant  layer.  We  call  the  scattered  wave  a  Stokes  wave  because  the 
scattering  process  is  a  three  wave  interaction  analogous  to  Raman  scattering.  The  Stokes 
wave  with  frequency  ids  =  lup  —  id  carries  the  information  encoded  on  the  signal  wave  back 
to  the  vehicle.  We  will  show  that,  whereas  much  of  the  scattered  Stokes  wave  propagates 
away  from  the  vehicle,  a  significant  fraction  is  returned  to  the  vehicle. 

What  is  remarkable  is  this.  The  ratio  of  the  power  flux  of  the  Stokes  wave  received 
at  the  vehicle  to  the  power  flux  contained  in  the  signal  wave  at  the  plasma  edge  can  be 
between  0.7  and  2  percent.  This  means  that  reception  of  GPS  signals  may  be  possible 
because  one  simply  needs  an  onboard  receiver  approximately  100  times  more  sensitive  than 
commercially  available  hand  held  receivers  or  use  sufficiently  larger  antenna.  We  shall  discuss 
in  the  Results  section  the  sensitivity  required  for  a  variety  of  sources. 

Communications  from  the  vehicle  requires  two  power  sources  on  the  vehicle.  One,  which 
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"Stokes”  wave 


Signal  wave 


Pump  wave 


FIG.  2:  The  concept  for  communication  from  the  vehicle.  Although  drawn  in  such  a  way  that 
the  angles  of  pumping,  Stokes  and  signal  waves  are  different,  the  optimal  configuration  is  when  all 
angles  are  the  same,  i.e.  Stokes  and  pump  waves  are  generated  in  the  same  direction  as  the  target 
of  the  desired  low  frequency  signal. 

we  term  the  Stokes  wave  generator,  will  also  carry  the  signal.  The  other  is  the  pump  wave. 
Both  have  carrier  frequencies  above  that  of  the  maximum  of  the  plasma  frequency.  Their 
nonlinear  interaction  in  the  plasma  produces  an  oscillations  of  frequency  u>  =  u>p  —  u>s- 
Consider  Figure  2.  For  zr  <  z  <  R  where  zT  is  determined  by  uji(zr)  =  u>  cos  <ft  and  <p  is 
calculated  from  the  differences  in  propagation  directions  of  the  pump  and  Stokes  waves,  the 
oscillation  does  not  propagate  and  its  strength  decays  away  from  the  vehicle.  Nevertheless 
this  oscillation  is  sufficiently  strong  to  act  as  a  power  source  for  a  propagating  wave  in  the 
region  z  <  zr  where  ui  cos  <fi  >  uL(z).  In  the  conclusion  we  analyze  what  power  is  required 
in  order  for  the  signal  to  be  detected  by  distant  receivers.  It  appears  that  even  if  we  use 
usual  available  on  the  market  generators  communication  can  be  put  into  practice. 


II.  OUR  RESULTS. 

A.  Equations  considered. 

We  begin  with  a  detailed  analysis  of  the  two  dimensional  propagation  and  interaction  of 
a  signal  wave  of  frequency  uj,  a  pump  wave  of  frequency  u>p  and  a  Stokes  wave  of  frequency 
ujs  through  a  plasma  with  a  given  density  profile  n0(z)  where  2  is  the  direction  normal  to 
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the  vehicle.  The  key  equation  is  a  modification  of  the  well  known  Ginzburg  equation  [4] 


d_ 

dz 


(  So  A  dH\  So  d^H_  Q2  -- 

\£(z,fl)J  dz  J  e(z,  f2)  dy2  c 2 


V  x 


(2) 


for  the  magnetic  field  amplitude  (H(y,z),  0, 0)e  lfit  of  an  oscillation  of  frequency  f2.  In  (2), 
the  effective  electric  susceptibility  is 


=  e0 


(,  u M 

\  W 


1 


1  +  i  v/Cl 


(3) 


(ujl(z)  is  the  local  plasma  frequency  and  v  the  collision  frequency).  The  susceptibility  is  due 
to  the  linear  response  of  the  plasma  to  the  electric  fields  of  whichever  waves  are  involved. 
The  nonlinear  current  j^L  will  be  determined  both  by  the  product  of  the  plasma  density 
distortion  with  the  linear  current  and  the  nonlinear  response  of  the  electric  velocity  field  due 
principally  to  dynamic  pressure  forces.  We  observe  that,  for  12  3>  maxwi(z),  the  electric 
susceptibility  is  approximately  £o  and  the  left  hand  side  of  the  nonlinear  Ginzburg  equation 
(2)  is  the  usual  wave  operator. 

How  do  we  use  (2)?  For  the  case  of  communication  to  the  vehicle,  we  use  it  in  two  ways. 
First  with  j^L  =  0,  we  determine  for  f2  =  u  and  H(y,z )  =  //(z)el^c^sin^,  the  field  H(z ) 
from  which  the  distortion  to  the  plasma  produced  by  the  incoming  wave  is  calculated.  In 
this  instance,  H(z)  satisfies 


d2H 


dz 2 


1  d£(z,w)dH 


£(z,w) 

ejz^j) 

£o 


+ 


dz  dz 
—  sin2  <p  ]  H  =  0. 


(4) 


A  glance  at  the  third  term  shows  that  propagation  is  impossible  for  e/eo  <  sin2  0  or,  from  (3), 
for  uj cos 0  <  ujl(z).  The  importance  of  the  resonance  layer  where  £(z,u>)  ~  0  is  seen  from 
the  denominator  in  the  second  term.  Having  solved  for  H(z)  from  (4)  we  can  then  calculate 
the  plasma  distortion  field  6n(z).  Its  interaction  with  the  pumping  wave  then  produces  a 
nonlinear  current  yv /.  which  gives  rise  to  the  Stokes  wave.  The  Stokes  wave  Hs{y,  z)  and  its 
propagation  is  calculated  by  solving  (2)  with  this  Jnl  and  appropriate  boundary  conditions 
at  the  plasma  edge  and  at  the  vehicle.  Our  goal  is  to  determine  Hs(y,z  =  R).  We  give 
the  results  of  both  the  numerical  simulation  and  an  analytic  estimation.  The  latter  takes 
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advantage  of  the  fact,  that,  for  the  Stokes  wave,  cus  maxc ol(z)  and  that  the  principal 

plasma  distortion  occurs  at  the  resonance  layer. 

For  communicating  from  the  vehicle,  we  solve  (4)  with  the  right  hand  side  given  by 
— V  x  ^fjNL  with  j/vL  calculated  from  the  nonlinear  interaction  of  the  pump  and  Stokes 
waves.  Here  the  goal  is  to  calculate  the  flux  of  power  of  the  signal  wave  with  frequency 
u>  =  Up  —  us  as  it  leaves  the  plasma  edge  in  the  direction  of  some  distant  receiver. 


B.  Analytic  estimations. 

1.  Model. 


We  shall  study  a  very  idealized  situation  when  the  plasma  sheath  is  a  flat  slab.  The 
plasma  density  is  a  linear  function  of  the  horizontal  coordinate  2 

2  +  L 


n0(z)  =  no 


R  +  L' 


(5) 


In  this  geometry  the  vehicle  is  the  vertical  wall  placed  at  z  =  R.  The  plasma  density  near 
the  vehicle  is  no-  The  plasma  contacts  the  vacuum  at  z  =  —L,  where  n  =  0.  We  shall 
study  two  situations:  communication  to  the  vehicle  and  communication  from  the  vehicle. 
In  both  cases,  three  almost  monochromatic  electromagnetic  waves  exist  in  the  plasma.  Two 
of  them  have  high  frequencies  up  (pumping  wave),  lus  (Stokes  wave).  The  third  one  has  low 
frequency  u,  satisfying  the  condition 


U)  =  UJP  -  US. 


(6) 


In  the  “to  the  vehicle”  case  u>  is  the  circular  frequency  of  the  incoming  signal.  In  the  “from 
the  vehicle”  case,  u>  is  the  circular  frequency  of  the  outgoing  signal.  In  both  these  cases,  the 
low-frequency  signal  plays  a  key  role.  Because  the  local  plasma  frequency  at  2  =  0  is  u, 


2  ^  no 

cj  = 


Mso  R  +  L 

Let  us  denote  also  the  Langmuir  frequency  at  the  vehicle  as 


(7) 


Wr  = 


e  »o 
Mso 


Thus 


R  +  L 


u> 


P 

fl 
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In  a  realistic  situation  fi  ~  9 GHz  (it  corresponds  to  no  =  1018m-3),  /  ~  2 GHz,  R  +  L  = 
1  m,  and  L  ~  0.05m.  The  wavelength  of  the  incoming  signal  in  the  vacuum  is  A  =  c/f  = 
0.15m,  so  that  A  >  L.  We  point  out  that  in  the  case  of  low-frequency  wave  reflection  from 
the  ionosphere,  the  situation  is  the  opposite  A  «  L. 

We  shall  assume  that  the  ions’  positions  are  fixed  and  the  plasma  is  cold  ( Te  ~  0).  The 
magnetic  field  has  only  one  component  Hx.  The  electric  field  has  two  components  Ey,  Ez. 
Neither  the  electric  nor  magnetic  fields  depend  on  the  ^-coordinate. 

2.  “To  the  vehicle.” 

We  would  like  to  estimate  the  ratio 

_  Ss(z  =  R) 

T'S  —  r, 

*->() 

of  the  fluxes  of  the  squared  scattered  field  to  the  squared  incoming  signal  field  and  express 
it  as  a  function  of  pump  power  flux  Sp  measured  in  Watts  per  square  meter. 

We  can  make  an  analytic  estimation  of  the  three-wave  process  efficiency.  The  main 
contribution  comes  from  the  vicinity  of  2  =  0.  The  reason  comes  from  the  fact  that  the  real 
part  of  dielectric  susceptibility  (3)  for  the  low  frequency  signal  wave  has  a  zero  at  this  point. 
It  means  that  our  nonlinear  current  on  the  right  hand  side  of  the  Ginzburg  equation  has  a 
very  sharp  peak  near  2  =  0.  A  typical  plot  of  right  hand  side  is  given  in  Figure  3.  This 


FIG.  3:  The  typical  right  hand  side  (absolute  value)  of  the  Ginzburg  equation  in  the  “to  the 
vehicle”  case.  Logarithmic  scale.  One  can  see  that  the  main  contribution  comes  from  the  region 
of  the  point  2  =  0. 
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issue  is  discussed  in  more  detail  in  Appendix  A  2. 

If  we  consider  a  high  frequency  pumping  wave  we  can  use  the  plane  wave  approximation 

Hp(y,z,t )  =  Hpei('~0Jpt+kpy~KpZK 

The  low  frequency  signal  wave  can  be  written 

H0(y,t)  =  H(y,z,t)  |2=0  =  H(z)\z=0ei{-“t+ky). 

For  the  Stokes  wave,  whose  frequency  is  higher  than  the  plasma  frequency,  one  can  use  the 
following  approximate  Ginzburg  equation 

Qz2  +  k2Hs  =  fs,  (8) 

where  fs  is  calculated  from  the  curl  of  the  nonlinear  current.  To  solve,  we  use  the  method 
of  variation  of  constants.  We  find 


Hs  =  Cxe'KsZ  +  C2e-'KSZ, 
C[eiKSZ  +  C'e~iKsZ  =  0, 

-L 


One  can  say  that  C\  is  the  amplitude  of  the  Stokes  wave  propagating  to  the  vehicle  and  C2  is 
the  amplitude  of  the  anti-Stokes  wave  propagating  from  the  vehicle.  The  main  contribution 
to  C\(R)  arises  from  the  vicinity  of  z  =  0,  where  fs(z)  is  almost  singular 

1  R 

Cl{R)  =  2i^/  fs(y)e~iKSVdy - 


i 

2i/c5 


+oo 

J  fs(y)e~lKsvdy. 

— oo 


After  some  simple  but  tedious  calculations  (see  Appendix  A  2)  one  finds 

Cx(R)  ~  27Ti-^  —  cos(2 e)  sin(0)//p//*(O).  (9) 

Me*  EqC 

where  9  is  the  pumping  incident  angle. 


9 


FIG.  4:  Dependence  of  Ci(f?)  on  the  signal  incidence  angle  <f>. 


Details  of  these  calculations  are  given  in  Appendix  A  2.  The  angular  dependence  of  H( 0), 
which  we  call  p(o)-  can  be  calculated  numerically  by  solving  the  homogeneous  Ginzburg 
equation.  In  Fig.  4,  we  plot  the  product  psincf)  against  (f>.  At  the  optimal  value  <t>  ~  0.5, 
p((j>)  sin  <p  ~  1/4, 


Using  the  expression  Sp  =  \Hp\2/(£qc),  one  gets 


Ms  = 
2 


Cl 


H 


C£  o 


4  (mc2)  £0C  COS  ^  lVUm-2 


\Wm~2 
S„ 


(10) 


(11) 


For  the  optimal  values  of  incidence  angles  (9  =  0,  0  ~  0.5),  the  given  plasma  parameters 
and  L  ~  0.05m,  one  gets  the  following  maximum  value  of  the  efficiency  coefficient 


ps  ~  0.9  x  10 


-li 


1  Wm~2 

This  is  consistent  with  what  we  obtain  by  direct  numerical  simulation. 


(12) 


3.  “ From  the  vehicle 


Equation  (2)  can  be  rewritten  in  the  following  form 

±idu  +  _  M'l h  =  JL  (Gmk) 

dze  dz  \£0  C2  £0 )  dz  \  £  ) 


£dy^jNL'z' 


(13) 
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It  is  not  too  surprising  that  that  the  dominant  contribution  to  the  RHS  of  (13)  is  the  first 
term  and  arises  from  the  neighborhood  of  z  —  0.  Again,  just  as  in  the  “to  the  vehicle”  case, 
the  resonant  layer  acts  as  a  transmitting  antenna  which  will  beam  the  message  contained  on 
the  Stokes  wave  to  a  distant  receiver  at  frequency  u>  =  uip  —  u>s.  the  dominant  contribution 
comes  from  the  first  term  on  the  RHS  of  (13)  and  from  the  neighborhood  of  z  =  0. 

If  we  consider  the  expression  for  (, ]NL)y  given  in  Appendix  A,  we  can  see  that  in  the 
case  cjo  a>s>  up  the  first  term  (A5)  is  the  major  one  in  the  vicinity  of  resonant  layer.  The 
resonant  layer  works  like  radiating  antenna. 

Using  the  simplified  nonlinear  current  expression  and  considering  the  pumping  and  Stokes 
waves  as  plane  waves  one  finds 


H(-L )  ~  -i 

M-L) 


etjgL  sin  0  1 
2M£qC?U>sOJp  A  p  s 


4>'\{—L )  +  iKo0i(— L) 


1  - 


g\Acos<})  _  j 

Acoscf) 


(14) 


Where  A  =  Lu)q/c. 

Using  the  solutions  of  the  approximate  homogeneous  equations,  we  can  estimate 
(t>'i(z)/<pi(z)\z=-L  —  1  /T.  Thus  for  k0T  =  A cos0  1,  one  finds 


H(-L)  ~ 


ecjpLsin^  1 

^MsocAusWp  A  p  s 


For  the  power  density,  we  have 

9=1 

32 


eL  \2  1 

M?j 


U— ^  sin2  4>SsSp. 


J 


(15) 


£oC  \ujsujP/ 

This  result  is  quite  clear  from  physical  point  of  view.  The  larger  <f>  is,  the  longer  is  the 
distance  over  which  the  signal  wave  is  generated  in  the  plasma. 

In  our  simulations,  A  ~  2.1  and  in  this  case  we  cannot  use  the  simplified  expression  given 
above.  Instead  we  find, 


x  tan2  0(1  —  2 


s  =  }l(iL)2- x 

8  \Mc2)  e0c  \ujsujpA/ 
sin(A  cos  0)  1  —  cos(A  cos  0) 


Acos0 


+  2 


A2  cos2  0 

1 


(16) 


-SsSp. 


1  +  CdeT  COS2  0~°  P 

Here  we  introduced  the  coefficient  Cder  =  («o0i/0i)2  the  value  of  which  we  obtain  from  our 
numerics. 
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Finally,  we  find 


(17) 


1  +  Oder  COS2  (j)^S^P' 


(18) 


The  subscripts  refer  to  the  frequencies  of  the  onboard  pump  waves.  Again,  we  find  the 
magnitude  and  angular  dependence  to  be  consistent  with  our  numerical  results. 

C.  Numerical  simulations. 

The  equation  we  solve  numerically  in  all  cases  is  the  Ginzburg  equation  (2)  including  all 
terms  on  its  right  hand  side.  The  boundary  conditions  are  given  aX  z  =  L\  =  —  L  —  (L  +  R), 
in  the  vacuum  beyond  the  plasma  edge  and  at  z  =  R,  the  vehicle. 

We  tested  the  robustness  of  the  code  by  allowing  for  both  finite  and  zero  conductivity 
of  the  vehicle  surface  in  the  “to  the  vehicle  case”.  During  the  simulation  in  the  “from  the 
vehicle”  case  we  also  redid  the  simulation  with  the  derivative  of  the  magnetic  field  at  the 
vehicle  equal  to  zero.  In  all  the  cases,  the  influences  of  the  differing  boundary  conditions 
were  negligible. 

1.  “To  the  vehicle.” 

As  the  first  step  in  the  “to  the  vehicle”  case  we  have  to  find  the  profile  of  the  incoming 
magnetic  field  in  the  plasma.  We  used  an  incident  angle  0  =  0.5.  It  will  be  shown  later 
that  this  angle  is  an  optimal  value  but  it  is  good  for  an  initial  evaluation  of  the  possibility 
of  communication.  We  consider  the  incoming  signal  as  a  monochromatic  plane  wave  of  a 
given  frequency  /0  =  2 GHz  and  amplitude  H0.  The  current  is  equal  to  zero.  In  this  case, 
the  boundary  conditions  are 


z 

FIG.  5:  Incoming  signal  magnetic  field  profile. 
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FIG.  6:  Incoming  signal  electric  field  (z-component)  profile. 


The  resulting  profile  of  the  magnetic  field  is  shown  in  Fig.  5.  The  profile  of  Ez(z )  is  shown 
in  Fig.  6.  At  the  next  stage,  we  consider  an  incident  low  frequency  magnetic  field  profile 
as  a  source  of  distortion  of  the  plasma  density  profile  and  take  into  account  currents  due  to 
the  presence  of  a  pump  wave.  The  pumping  wave  angle  6  =  0.0.  Our  goal  is  to  calculate  the 
scattered  field  Hs  with  frequency  us  =  —  u.  In  this  case,  the  boundary  conditions  are 


z 

z 


=  -Lu  ^  +  iKSH  =  0, 


(21) 

(22) 


The  profiles  of  the  magnetic  fields  Hs  for  two  different  pumping  frequencies  are  shown  in 
Figures  7  and  8.  We  note  that  the  resonant  layer  2  =  0  acts  as  if  it  were  a  source.  In  the 
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FIG.  8:  Magnetic  field  profile  of  the  Stokes  wave.  Pumping  frequency  18GHz. 

“to  the  vehicle”  case,  it  is  convenient  to  introduce  the  function  ns  as  the  ratio 

Ss{z  =  R) 


Ms  = 


So 


of  the  scattered  field  flux  to  the  incoming  signal  flux  and  express  it  as  a  function  of  pump 
flux  Sp  measured  in  Watts  per  square  meter.  We  found 


5 

up  =  2ir  *  12 GHz,  max(fis)  —  2.2  x  10~12  p — 

l\Vm~z 

Q 

Up  =  2n  *  18 GHz,  max(/is)  —  0.63  x  10~n  p  9. 

*  1  Wm~2 

These  results  are  in  a  good  agreement  with  the  analytic  estimation  (12).  Any  difference  is 

due  to  the  fact  that  the  pumping  frequency  is  not  sufficiently  high  to  neglect  the  plasma 
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FIG.  9:  Dependence  of  power  conversion  efficiency  coefficient  fis  on  angles.  “To  the  vehicle”. 
Pumping  frequency  12GHz. 


1.8 


FIG.  10:  Dependence  of  power  conversion  efficiency  coefficient  fis  on  angles.  “To  the  vehicle”. 
Pumping  frequency  18GHz. 

frequency.  The  reason  we  used  these  frequencies  and  not  much  higher  ones  was  that  they 
are  available  on  standard  microwave  equipment  and  devices. 

In  order  to  investigate  the  dependence  of  the  result  on  the  angles  0,  9P ,  we  calculated 
//  for  various  different  choices.  The  results  are  shown  in  Figs.  9-10. 

As  one  can  see,  in  the  “to  the  vehicle”  case  we  have  a  very  good  agreement  between 
the  analytically  estimated  angular  dependence  (11)  and  the  numerical  results.  Namely,  we 
have  a  maximum  at  pumping  angles  close  to  9  =  0  and  the  efficiency  coefficient  jis  goes 
to  zero  at  the  vicinity  of  9  =  7r/4  in  a  agreement  with  the  cos(20)  dependence.  So  we  can 
formulate  a  simple  rule:  in  order  to  get  the  best  possible  performance,  send  the  pump  wave 
in  a  direction  perpendicular  to  the  plasma  edge  surface. 
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FIG.  11:  Dependence  of  the  power  conversion  efficiency  coefficient  /j,s  on  several  pumping  angles, 
in  the  “to  the  vehicle”  case.  Pumping  frequency  12GHz. 


FIG.  12:  Dependence  the  of  power  conversion  efficiency  coefficient  /<s  on  several  pumping  angles, 
in  the  “to  the  vehicle”  case.  Pumping  frequency  18GHz. 

2.  “From  the  vehicle.” 

In  the  “from  the  vehicle”  case  we  calculate  the  magnetic  field  of  the  low  frequency  wave 
generated  by  plane  pump  and  Stokes  waves.  Following  the  optimal  strategy  in  this  case, 
described  in  the  analytic  part  of  the  paper,  we  take  all  angles  equal  to  each  other  (ft  =  9  = 
7r/4.  In  this  case,  the  boundary  conditions  are 


dH  . 

2  =  -L 1,  —  +  IKqH  =  0, 
oz 

(23) 

z  =  R,  H  =  0. 

(24) 

Here  H(z)  is  the  magnetic  field  of  the  signal  wave  with  frequency  u  =  up— us-  The  boundary 
condition  at  z  =  R,  H  =  0,  gives  us  the  worst  of  all  cases  by  definition. 
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FIG.  13:  Generated  low  frequency  magnetic  field.  Pumping  frequency  12GHz. 


FIG.  14:  Generated  low  frequency  magnetic  field.  Pumping  frequency  18GHz. 

The  low  frequency  magnetic  fields  for  two  different  pumping  frequencies  are  shown  in 
Figs.  13  and  14.  In  the  “from  the  vehicle”  case,  we  calculate  the  ratio 

_  S0ut(z  ~  L) 

M  “  SsSp 

of  the  output  signal  flux  to  the  product  of  the  pump  and  Stokes  fluxes  and  express  it  as  a 
function  of  the  optimal  angle. 

We  found 

up  =  2n  *  12 GHz,  maxf/j)  ~  1.8  x  10-16— 1 — 

1  Wm~l 

up  =  2n  *  18 GHz,  ma x(/z)  ~  3.0  x  10”17-— - — -. 

1  Wm~z 

In  order  to  investigate  the  dependence  of  the  result  on  the  angles  0,  9P ,  Os,  we  calculated 
/j  for  various  different  choices.  The  results  are  shown  in  Figs.  15-16. 
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FIG.  15:  Dependence  of  the  power  conversion  efficiency  coefficient  //  on  optimal  angle.  “From  the 
vehicle”.  Pumping  frequency  12GHz. 


FIG.  16:  Dependence  of  power  conversion  efficiency  coefficient  fi  on  optimal  angle.  “From  the 
vehicle” .  Pumping  frequency  18GHz. 

In  the  “from  the  vehicle”  case,  the  situation  is  even  simpler.  As  it  was  shown  in  Section 
II B  3  the  power  conversion  is  optimal  if  we  radiate  both  the  pump  and  Stokes  waves  in  the 
direction  of  the  desired  signal  wave  propagation.  The  estimated  angular  dependence  (16) 
can  be  fitted  with  good  accuracy  to  the  numerical  results  using  only  one  tuning  coefficient 
Cder ■  It  is  shown  that  this  coefficient  weakly  depends  on  the  pumping  frequency. 


III.  PRACTICAL  USAGE  OF  THE  OBTAINED  RESULTS. 

Let  us  now  discuss  the  practical  usage  of  this  approach  for  receiving  at  and  transmitting 
from  the  vehicle.  For  the  “to  the  vehicle”  case  we  consider  the  problem  of  receiving  even 
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GPS  signals.  Let  us  estimate  the  resulting  attenuation  coefficient.  Given  a  pump  waveguide 
aperture  of  3 cm  x  3 cm  and  a  pump  power  of  3 kW,  this  gives  Sp  =  3.3  x  10 6Wm~2.  One 
can  use  the  pulse  regime.  In  this  case,  even  for  pulses  10-3s  long,  every  pulse  still  contains 
more  than  106  periods  of  the  low  frequency  signal  and  we  can  get  much  higher  power  flux 
gpuUe  _  3  3  x  io 9Wm~2.  It  gives  us  the  attenuation  coefficients  nsSpulse 

HSS%ulse  ^  0.73  x  10“2,  up  =  2tt  *  12 GHz. 

HsS^lse  ~  2.1  x  10“2,  up  =  2n  *  18 GHz. 

The  usual  level  of  a  GPS  signal  at  the  Earth  surface  is  about  —127.5 dBm  (1  Decibel  per 
milliwatt  is  equal  to  1  dBm  =  10  log10(P/lmtP)).  Indoors,  one  must  use  high  sensitivity  GPS 
receivers.  Many  general  purpose  chipsets  have  been  available  for  several  years.  Presently, 
the  market  offers  sensitivities  —157.5 dBm  (for  example  [5]).  Using  the  definition  of  dBm 
one  can  see,  that  it  is  possible  to  receive  a  signal  with  an  attenuation  about  10-3.  Also  it  is 
possible  to  use  a  much  bigger  antenna  on  the  vehicle  than  in  the  case  of  a  handheld  device. 
In  this  case,  it  is  even  possible  to  receive  a  signal  using  the  continuous  rather  than  pulsed 
regime  for  a  klystron  pump.  So  even  at  the  angles  far  from  optimal,  one  can  receive  GPS 
signals.  Further,  we  used  maximum  value  of  the  plasma  thickness.  If  the  plasma  sheath  is 
thinner,  the  angular  dependence  is  broader. 

Some  characteristics  of  klystron  amplifiers  available  on  the  open  market  are  given  here 
[5,61 _ 


Model 

Frequency  (GHz) 

Power  (kWatt) 

Mass  (kg) 

LD4595 

14.0-14.5 

3 

40 

LD7126 

17.3-18.4 

2 

27 

In  the  “from  the  vehicle”  case,  because  of  sensitive  land  based  receivers,  all  we  need  is  to 
have  a  reasonable  signal.  Let  us  estimate  an  incoming  power  on  some  land  based  antenna. 
First  of  all,  for  any  real  antenna  we  have  to  take  into  account  the  decrease  of  a  signal  due 
to  diffraction  broadening.  If  the  diameter  of  the  land-based  antenna  (Figure  17)  is  D0 ,  the 
diameter  of  the  signal  flux  after  some  long  distance  l  will  be 

D(D  -  (25) 

It  means  that  if  we  have  power  flux  at  an  antenna  Sa,  the  power  flux  at  the  edge  of  the 
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FIG.  17:  Schematic  plot  of  beam  diffraction. 


plasma  after  a  distance  /  will  be 

*-*(£)’•  <26) 

For  example,  for  an  antenna  of  the  diameter  equal  to  5m,  after  100km 

S0  ^  1.1  x  10_55^. 

Now  one  can  calculate  the  sensitivity  of  the  receiver  needed.  Let  us  suppose  that  the 
signal  beam  outgoing  from  the  vehicle  has  diameter  D0  =  lm,  signal  frequency  /  =  2 GHz 
and  corresponding  wave  length  A  =  1.5  x  10_1m,  the  land  based  antenna  has  a  diameter 
Dlb  =  5 m  and  is  situated  at  a  distance  l  =  lOOfcm.  Using  the  previous  results  for  diffraction, 
the  pumping  klystrons’  powers  from  the  table  above  and  the  expression  Sout  =  I^SpSs,  one 
can  get  for  the  power  on  the  land  based  receiver 

2 Diy 

l\  ) 

We  now  list  for  two  different  frequencies,  the  corresponding  powers  in  Watts  at  the  receiving 
antenna. 


1.8  x  10-8S, 


out • 


(27) 


up  =  27t  *  12 GHz, 

PA  ~  1.8  x  10'8  *  1.8  x  1(T16  *  9  x  106VUm-2  *  25m2 
~  0.73  x  10-15W; 
up  =  2n  *  18 GHz, 

PA  ~  1.8  x  10-8  *  3.0  x  10-17  *  4  x  1012VUm-2  *  25m2 
~  0.54  x  10_17W. 
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The  GPS  receiver  mentioned  above  has  a  sensitivity  about  —  160d£?m  ~  10_19IV.  Even 
with  such  a  modest  size  of  the  antenna  and  ordinary  klystrons  one  can  receive  the  signal  at 
almost  any  angle. 

As  a  final  remark  one  can  conclude  that  proposed  method  for  communication  with  and 
from  the  supersonic  vehicle  is  realistic  even  using  standard  devices  available  on  the  open 
market. 

IV.  CONCLUSION. 

We  overcame  the  challenges  described  in  the  begining  of  this  report.  We  achieved  the 
complete  success  in  estimation  of  applicability  of  the  method  we  have  proposed.  Current 
technique  can  be  utilized  for  communication  with  hypersonic  vehicles  even  using  current 
microwave  devices  available  on  the  open  market. 

Following  “Report  of  Department  of  Defense,  Advisory  Group  on  Electron  Devices.  Spe¬ 
cial  Technology  Area  Review  on  Vacuum  Electronics  Technology  for  RF  applications”  at 
least  since  year  2000  much  more  powerfull  and  advanced  klystrons  are  available  for  military 
purposes.  It  makes  our  estimations  even  more  optimistic. 

The  code  we  created  can  be  easyly  used  for  simulations  with  an  arbitary  plasma  density 
profile. 

Results  of  the  current  work  were  published  [7]  in  the  Journal  of  Apllied  Physics,  which 
is  a  journal  of  such  a  respectable  organization  as  American  Institute  of  Physics. 

APPENDIX  A:  RIGHT  HAND  SIDE  OF  THE  GINZBURG  EQUATION. 

1.  General  case. 

Consider  the  Ginzburg  equation  for  a  wave 

H3(y,z,t)  =  H3(z)e~i“3t+ik3y, 

and  calculate  the  right  hand  side  of  (2)  in  terms  if  the  fields  Hi(y,  z,  t ),  H2(y,  z ,  t).  In  the  “to 
the  vehicle”  case,  H i  will  represent  the  pump  wave,  Hi  the  signal  wave  and  H3  the  Stokes 
wave.  In  the  “from  the  vehicle”  case,  H3  will  be  the  signal  and  H\  and  H2  the  pump  and 
signal  carrying  Stokes  waves  respectively.  In  all  cases  u>3  =  u2  —  u>i ,  k3  =  k2  —  k\ .  We  find 
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2.  Approximate  right  hand  side.  “To  the  vehicle”  case. 


In  the  “to  the  vehicle”  case,  the  main  contribution  comes  from  the  terms  containing  poles 

[V  x  jNL(ff,Hp,k,kp,ks,(J,(jp,us)Jx  ^  (A7) 

_ ±(i 1A^.)h*H  + 

M2  (1  -  iv/u)p)  u3ep  (1  +  w/ujp)  UJ2  dz  \  £*  )  p 

_ 1  dHp  k  d  /  d_  ( n0(z)\\  _ 

M2  uj2£p  (1  +  w/ijjp)  dz  u3  (l —w/u>)  dz  \  dz  y  £*  )) 

Moe3  k  d  ( n0{z)\ 

M2  u)3  (1  -  w/ui)  dz  y  e*  )  p 
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Assume  that  the  high  frequency  pumping  wave  remains  undisturbed.  Then 

Hp(y,z,t)  =  Hpel('kpy~KpZ~wt\ 

we  find 

fs(z)  =  -[VxjNL]xe-^z 

and 

1  R 

Ci(R)  =  j —  /  fs(y)e-iKsydy  * 

— L 

+oo 

— oo 

After  several  integrations  by  parts  in  the  second  term  of  [v  x  jNL ]  ,  taking  into  account 
ks  =  kP  —  k,  one  finds 


—ik  p3H  H* 

Ci(R)  —  — - ~  .9.  .s'  ( kp(kp  —  £5)  +  (kp+ 


k5)k, 


2ks  M2£0uj2lj3 

'2X 


(A9) 


_  ^  7  —  (' 

c2  J  dz{ 


dz. 


For  most  pumping  angles,  and  using  the  fact  that  cop  one  can  substitute  cjg  ~  up 

and  consider  incidence  angles  of  pumping  and  Stokes  wave  to  be  close  in  absolute  value. 
Following  Fig.  1  the  pumping  incidence  angle  is  0  and  low  frequency  signal  incidence  angle 
is  4 >• 


Ci(R) 


— ik  e3Hp  up 
2 ns  M2£0ujIuj3  c2 


cos(20) 


(A10) 


\  £*  J 

Using  integration  by  parts  once  more  one  can  get 

(kp  +  k5)  e3Hp 


Ci(R)  ~ 


2ks  M2£qC3OJ2 

x  T:mz) 


cos(20)  sin(0)  x 


H*e-'l{KS+^)zdz. 


(All) 


Calculating  this  integral  by  residues  and  taking  into  account 

.  .  .  .  L  J1 

n0(z)  ~  n0{z  =  0)  =  n0r  =  n0  —  , 

L  +  it 

£0  _  L  _  .  v 

.  -r)  d  i  L, 
£  z  + 10  u> 
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we  finally  get 


Cx(R) 


2 "  .X w  cos<28)  ™WW{0)  = 

0  L 

=  27ri  —  cos(20)  sin(4>)HpH*(0). 

Me2  e0c 
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